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ABSTRACT 

While the use of galaxy clusters as tools to probe cosmology is established, 
their conventional description still relies on the spherical and/or isothermal mod- 
els that were proposed more than 20 years ago. We present, instead, a depro- 
jection method to extract their intrinsic properties from X-ray and Sunyaev- 
Zel'dovich effect observations in order to improve our understanding of cluster 
physics. First we develop a theoretical model for the intra-cluster gas in hydro- 
static equilibrium in a triaxial dark matter halo with a constant axis ratio. In 
this theoretical model, the gas density profiles are expressed in terms of the in- 
trinsic properties of the dark matter halos. Then, we incorporate the projection 
effect into the gas profiles, and show that the gas surface brightness profiles are 
expressed in terms of the eccentricities and the orientation angles of the dark 
halos. For the practical purpose of our theoretical model, we provide several em- 
pirical fitting formulae for the gas density and temperature profiles, and also for 
the surface brightness profiles relevant to X-ray and Sunyaev-Zel'dovich effect 
observations. Finally, we construct a numerical algorithm to determine the halo 
eccentricities and orientation angles using our model, and demonstrate that it 
is possible in principle to reconstruct the 3D structures of the dark halos from 
the X-ray and/or Sunyaev-Zel'dovich effect cluster data alone without requiring 
priors such as weak lensing informations and without relying on such restrictive 
assumptions as the halo axial symmetry about the line-of-sight. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general 
— X-rays: galaxies: clusters 
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1. INTRODUCTION 

Understanding of statistical properties of dark matter halos has been significantly ad- 
vanced in recent years, largely owing to the development of high-resolution numerical simu- 
lations (Navarro, Frenk, & White 1996, 1997; Fukushige & Makino 1997; Moore et al. 1998; 
Jing & Suto 2000, 2002). Given those theoretical/empirical successes, a next natural ques- 
tion is how to apply them for the description of real galaxy clusters. While there exist a 
number of attempts along this line, they usually make the unrealistic assumption of the halo 
spherical symmetry (Makino, Sasaki & Suto 1998; Suto, Sasaki & Makino 1998; Yoshikawa 
& Suto 1999). This paper describes a methodology to reconstruct the three-dimensional (3D) 
structure of dark halos from the two-dimensional (2D) surface brightness profiles of intra- 
cluster gas from the X-ray and/or Sunyaev-Zel'dovich (SZ) effect observations assuming that 
the halos are triaxial ellipsoids with a constant axis ratio. 

Indeed it is a classical problem in astrophysics to determine the 3D properties of astro- 
nomical objects from the observed 2D counterparts (e.g., Lucy 1974; Fabricant et al. 1984; 
Dehnen & Gerhard 1993; Binney, Davies, & Illingworth 1990; Gerhard & Binney 1996). De- 
projecting galaxy clusters is one of the most important applications. Fabricant et al. (1984) 
analyzed the X-ray surface brightness map of clusters, and showed that their mass distribu- 
tion is far from being spherical. Zaroubi et al. (1998) developed a general method of depro- 
jecting the 2D images of rich clusters based on the Fourier slice theorem to reconstruct the 
3D cluster structures. Later their technique was tested against numerically simulated galaxy 
clusters (Zaroubi et al. 2001). Yoshikawa & Suto (1999) proposed a deprojection method 
for spherical clusters based on the Abel integral. Reblinsky (2000) provided a parameter- 
free Richardson- Lucy algorithm to reconstruct the 3D halo potential, and demonstrated its 
stability by applying it to gas-dynamical simulations. Dore et al. (2001) provided a pertur- 
bative approach to the cluster deprojection, taking into account the non-isothermality and 
asphericity of galaxy clusters. Recently, Fox & Pen (2002) also considered the problem of 
deprojecting aspherical clusters. They first constructed a parameterized 3D axisymmetric 
cluster model, and determined the 3D cluster shapes through the x 2 -fitting between the 
model predictions and the simulated cluster data. 

All the previous approaches, however, were based on rather restrictive assumptions 
such as the cluster axial symmetry (e.g., oblate or prolate) about the line-of-sight, and/or 
the isothermality of galaxy clusters. Furthermore in their approaches, it was required to 
combine X-ray and/or SZ effect data with weak lensing (WL) map, which significantly limits 
the applicability of the previous methods. In fact, it is generally believed that it would be 
a difficult task to deproject the clusters from the 2D projected X-ray or SZ observables 
alone without having such restrictive assumptions and priors given the degeneracy of the 
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parameters due to the the projection process itself. 

Lee & Suto (2003; hereafter Paper I) suggested that one may understand the properties 
of the dark matter halos from the intra-cluster gas profiles. Assuming that the intra-cluster 
gas is in hydrostatic equilibrium in the triaxial dark matter halos, Paper I derived the 
3D density and temperature profiles of the intra-cluster gas from the 1st principles using 
the perturbation theory, and found an analytic relation between the eccentricities of the 
intra-cluster gas and the underlying dark halo (see eq.[28] in the Paper I). However, the 
perturbation result of the Paper I is valid only in the asymptotic limit of the low asphericity, 
and did not include the effect of the projection of the gas profiles on the plane of the sky. 
Here, we generalize the works of the Paper I to the case of highly aspherical clusters, and 
attempt to find a general relation between the dark halos and the intra-cluster gas taking 
into account the parameter-degeneracy caused by the projection. For the practical purpose 
of our theoretical modeling, we provide a series of empirical fitting formulae for the 2D and 
3D gas profiles which can be used as templates in comparing with the observed profiles of 
clusters, and demonstrate the degree of the feasibility of the halo reconstruction using the 
template formulae and the simulated surface brightness maps of X-ray and SZ effect. 

The organization of this paper is as follows. In §2 we construct a theoretical model 
for the intra-cluster gas in hydrostatic equilibrium within the gravitational potential of a 
triaxial dark matter halo with a constant axis ratio, and provide a series of empirical fitting 
formulae for the 3D gas profiles in a systematic manner. In §3 we consider the projection 
of the 3D gas profiles onto the plane of the sky, and provide another series of empirical 
fitting formulae for the 2D surface brightness density profiles. In §4 we describe a numerical 
algorithm to determine the 3D structures of the dark halos from the observed 2D cluster 
surface brightness profiles, and test the algorithm against a numerical toy model. Finally we 
discuss the final results, and draw our conclusions in §5. 

2. MODELING INTRA-CLUSTER GAS PROFILES 

2.1. Gravitational Potential of Dark Matter Halos 

To predict theoretically the X-ray and SZ profiles of galaxy clusters, one first needs a 
good physical model for the intra-cluster gas. Here we assume that the intra-cluster gas 
is either isothermal or polytropic, and in hydrostatic equilibrium within the gravitational 
potential generated by a concentric and coaxial triaxial dark matter halo. 

Consider a triaxial dark halo whose iso-density surfaces are given by the following equa- 
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tion: 



l-el l-e c 2 ' 

where the Cartesian system of coordinates (x, y, z) is aligned with the halo principal axes, 
oriented such that the rr-axis and z-axis run along the major and minor principal axes, 
respectively. The major axis length of the iso-density surface is denoted by R, while e& and 
e c (e?, < e c ) represent the two constant eccentricities of the ellipsoidal dark halos. Equation 
(1) implies that the density profile of a triaxial dark matter halo should be a function of the 
major axis length R. 



p(R) = ^,^r crit ^^ (2) 



We adopt the density profile of a triaxial dark halo proposed by Jing & Suto (2002): 

^cPcrit 

(R/Ro) a (1 + R/Ro) 

where Ro is the scale length, 5 C is the dimensionless characteristic density contrast with 
respect to the critical density p cr \t of the universe at the present epoch, and a represents the 
inner slope of the density profile. Jing & Suto (2002) showed that asilon the cluster scale 
and a 3/2 on the galaxy scale. For simplicity and definiteness, we focus on the case of 
a = 1 throughout this paper. 

The gravitational potential due to the ellipsoidal halos described by equation (2) is 
formally expressed as (Binney & Tremaine 1987): 



*( r ) = -^d-^d-^T +i y°°»-y' dr, (3) 

V Jo V( T + 1 )( r + 1 - e b)( r + 1 - e c) 

r. )ji 2 2 2 

«m) = 2 1 P (RWR, m * = JL- i + -^ + -L_. (4) 

Equations (3) and (4) along with equation (2) allow one to compute the triaxial halo grav- 
itational potential at least numerically. The halo gravitational potential $(r) depends on 
5 C , p CT it,-Ro as we ll as e b, e c , a. However, the dependence of $(r) on 5 C , p crit , and R can be 
separated out by introducing a dimensionless potential $ that depends only on eb, e c and a 
such that 

\ $(r;e 6 ,e c ,o;,5 c ,p cr it,-Ro) ,_v 



2.2. Gas Density and Temperature Profiles in terms of Halo Potential 

To determine the intra-cluster gas profiles in terms of the halo potential, it is necessary 
to specify the equation of state for the intra-cluster gas. We consider both the isothermal 
and the polytropic cases in order. 
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2.2.1. Isothermal gas 
The equation of state for the isothermal gas is given as 

P g (r) = P g0 - 



PgO 



(6) 



where P g and p g represent the gas pressure and the gas density, respectively. In what follows, 
the subscript of some physical variable indicates the value of that physical variable at the 
center r = 0. 



The density profile of the isothermal gas in hydrostatic equilibrium is given as 

Pg(r) = p g0 exp -«{4>(r) - $ } 
We define a dimensionless isothermal gas constant k as 



K 



A7iGfi g m p 5 c p cri tRl 
k B T g 



(7) 



(8) 



where m p is the proton mass, \i g is the mean molecular weight of the intra-cluster gas, k B is 
the Boltzmann constant, and T g is the (constant) gas temperature. Introducing P$(r): 



Pj, (r) = exp 



-{$(r)-$ } 



one can rewrite equation (7) as 



p g (r)=p g0 [F*(r)]« 



(9) 



(10) 



2.2.2. Polytropic gas 
The equation of state for the polytropic gas with the polytropic index of 7 (7^ 1) is given 



as 

and its density and temperature profiles in hydrostatic equilibrium are found as 1 

1/(7-1) 



(11) 



p g (r) = p g0 1 - « P {^(r) - $0} 



P g — P g o 



1 - Kp {$(r) - $0} . (12) 



1 Note that the definition of <& in Paper I is different from that given here by a constant offset. 
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We define a dimensionless polytropic gas constant 

_ 7 - 1 4:irGn g m p 5 c p crit Rl 



(13) 




(14) 



P g (r) = p g0 [1 - Kp + K p F$(r)} 



1/(7-1) 



T g (r) =T g0 [l-Kp + KpF*(r)]. 



(15) 



2.3. Empirical Fitting Formulae for the Gravitational Potential 

We have shown in §2.2 that the potential function F$(r) defines the density and the 
temperature profiles of intra-cluster gas completely. In general, F$(t) does not have a closed 
analytic form when the dark matter halos are triaxial ellipsoids. In Paper I, we computed 
F$(r) analytically with the perturbation theory assuming ej < « 1, and related the 
iso-density surfaces of the intra-cluster gas to those of the dark halos; if the halo iso-density 
surfaces are triaxial ellipsoids with the eccentricities of e a (a = b,c), then the iso-density 
surfaces of the intra-cluster gas is also well approximated as triaxial ellipsoids with the 
eccentricities of e a that are related to e a by 



where u = |r|/i? - Note that the right-hand-side of equation (16) are written in terms 
of the rescaled spherical radius u instead of the rescaled major axis length, which results 
from the fact that we neglected the higher-order terms 0(e£) in deriving equation (16) with 
the perturbation theory. Although equation (16) seems complicated, it is a slowly varying 
function, and well approximated as a constant in the range of < u < 1 (see Fig. 3 of Paper 
I) such that 



Strictly speaking, equation (17) is valid only in the limit of < C 1 over the 
range of < u < 1. To proceed further in a tractable fashion, however, we extend the 
validity of equation (17) to e 2 < e 2 < 1 over the whole range of u. In other words, we 
assume that the iso-density surfaces of the intra-cluster gas are triaxial ellipsoids with the 
constant eccentricities of e&, e c that are related to the halo eccentricities e^, e c by equation 
(17). Actually this turns out to be a good approximation as will be shown below. 



e 2 a _ 6(1 + u) ln(l + u) + u 3 - 3u 2 - 6u 
e| ~ 2u 2 [(l + u) ]n(l + u)-u] 



(16) 



% » 0.7 2 . 

4 



(17) 
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This assumption allows us to write the potential function F$(r) in terms of the major 
axis length of the gas iso-density surfaces, £, defined as 



Ux + 1 - e 2 + 1 - e c 2 i? 2 ^ + 1 - e 2 b + 1 - r» ' ' 



After some trials and errors, we find that the following empirical formula fits both equations 
(9) and (14) quite well: 

1 + rjf * \ q 



F *®=[TTW')' (19) 

where (3, p, q, and 77 are free parameters. The four free parameters are not constant but 
supposed to be functions of the halo eccentricities. We determine the functional form of each 
parameter by taking the following steps. 

1. Compute F$(t) numerically using equations (3) and (4) for various values of e& and e c . 

2. Fit those numerical data points to the empirical model (eq. [19]) and obtain the best-fit 
values of /3, p, q, and rj using the Levenberg-Marquardt method (Press et al. 1992). 

3. Finally model those best-fit parameters as functions of and e c . In practice, we find 
that all the parameters can be written as functions of a single variable fi = ef + e^. 

We find the following polynomials to provide good fits: 

(3 = c^o + c^i/i, (20) 

p = const., (21) 

q = c q0 + Cgi/i + c g2 /i 2 , (22) 

r] = c v0 + Cnifi. (23) 

Table 1 lists the best-fit values of the polynomial coefficients. The best-fit values of the 
constant p is determined to be unity for all cases. Figure 1 plots (3, q, and n as a function 
of pi. The filled circles indicate the best-fit values of the parameters, while the solid lines 
represent their polynomial fits. Note that 77 is one order-of-magnitude smaller than f3 for all 
cases, suggesting that the term associated with 77 can be safely neglected except at the outer 
part of clusters. The accuracy of the fits is illustrated in Figure 2 where is plotted 

against £. The solid lines represent the numerical results while the dashed lines correspond 
to equation (19) with the best-fit coefficients of the polynomials plotted in Fig. 1 and listed 
in Table 1. Our empirical fitting formulae reproduce the numerical results very well for all 
cases over a wide range of £ within a fractional error less than 0.5%. 
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The characteristic feature of our fitting model is that all the parameters are expressed 
in terms of fj, = ef + e\. Given that the first-order perturbation approach of the Paper I 
indicated the dependence of the halo potential on e\ + e 2 , one may have expected the fitting 
model to depend on e\ + e 2 rather than ef + e\. As a matter of fact, we indeed attempted 
first to model the parameters in terms of e\ + e 2 . But, it turned out that the parameter 
values do not scale in terms of e\ + e 2 c while they show very good scaling feature in terms 
of e\ + e\ as shown in Figure 1. We have not yet completely understood the origin of this 
scaling, but the result is empirically quite robust. 

In addition, we note that in the Paper I we derived the gas-halo eccentricity relation 
(eq. [16]) by using the first-order perturbation theory, and showed that the approximation 
error can be expressed as a function of e\ and e\ (see eqs [29] and [30] in the Paper I). The 
dependence of the approximation error on e\ and e\ implies that the higher order terms 
neglected in the first-order perturbation result may be scaled in terms of e\ and e\. Now 
that we have used equation (16) to construct the fitting model (eq. [19]), one may expect 
that equation (19) has similar scaling. In other words, we suspect that the dependence of 
equation (19) on e\ + e\ rather than e\ + e 2 may be related to the scaling of the higher order 
terms neglected in the first-order perturbation result of the Paper I. 

Note also that in our model the gas density/temperature profiles do not approach zero 
even at large u. This feature is ascribed to the hydrostatic equilibrium condition itself. In 
reality, the hydrostatic equilibrium condition may be satisfied only within some radius, e.g., 
the virial radius of the cluster. In fact, the halo density profiles at those regions are not 
well approximated by equation (2) because of the presence of substructures, and so on. This 
problem, however, is unlikely to limit the validity of our methodology in practice since the 
X-ray and SZ fluxes from regions beyond the cluster virial radius are usually negligible. 



Now we are in a position to model the X-ray and SZ surface brightness profiles of triaxial 
galaxy clusters by incorporating the projection effect into the 3D models developed in §2. 
Just for comparison, let us recall that the surface brightness profiles of spherical clusters can 
be readily evaluated as 



where 6 is the angular radius from the center of the cluster, d A is the angular diameter 
distance to the cluster, and L is the emissivity given in terms of p g and T g such that L oc 



3. 



PROJECTED SURFACE BRIGHTNESS PROFILES OF 
INTRA-CLUSTER GAS 




(24) 
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p 2 T g ' for bolometric X-ray and L oc p g T g for SZ observations, respectively. For the case 
of a triaxial cluster, however, the result becomes much more complicated because the 2D 
projection depends on the relative direction of the line-of-sight of an observer with respect 
to the principal axes of a triaxial halo. The projection effect of triaxial bodies on the plane 
of the sky is fully discussed by Stark (1977) and Binney (1985). In what follows we adopt 
the notation of Binney (1985). 

Let (9, <p) be the polar angle of the line-of-sight in the halo principal coordinate system 
(x,y,z), and let the observer's coordinate system be defined by Cartesian axes (x',y',z') 
with z'-axis aligned with the line-of-sight direction and rr'-axis lying in the (x, y) plane (see 
Fig. 1 in Oguri, Lee, & Suto 2003). Then the halo principal coordinate system (x,y,z) is 
related to the observer's coordinate system (x',y',z') by 

x = — sin <f)x' — cos (f) cos 9y' + cos (f) sin 9 z', (25) 
y = cos(f)x' — sin (f) cos 9y' + sin <p sin 9 z' , (26) 
z = sin9y' + cos9z. (27) 



If we use the major axis length £ defined in equation (18) to characterize the iso-density 
surfaces, then the projection of L onto the plane of the sky is written as 

/OO cy poo 

L{e)dz' = -= / L(z" 2 + \ 2 )dz", (28) 
-oo V J JO 



where 



z" = V / 7^' + ^), (29) 



A 



-L(Ax' 2 + Bx'y' + Cyy/ 2 , (30) 



sm 2 \ cos 2 9 

i. 



f = s in 2 9^cos 2 ( f>+ T -j I J+ T —t, (31) 
g = sin9sin2<p( -A x' + sin29 ( -cos 2 <p- ^^)y', (32) 



i cos 2 # / . 2 cos 2 <A sin 2 # 



B = cos 6 sin 20 1 A -, (34) 
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Hence T,(x',y') = X(A). In other words, the 2D isophotal curves of triaxial galaxy clusters 
are concentric and coaxial ellipses if their eccentricities e a of the gas iso-density surfaces 
are constants. As we noted before, however, it is not exactly the case for our model where 
the eccentricities e a of the halo iso-density surfaces are constant. Nevertheless this holds 
approximately (see eq.[17]). 

Let us define a dimensionless surface brightness E = S/S where S = S(0) and 
A = \/Rq. In the same spirit as for the 3D model (19), we propose the following empirical 
model for the 2D surface brightness profiles: 

S(A; e 6 , e c , a;, 7) = ( * + ^~ ) . (36) 

V 1 + CA s y 

Note that S depends not only on \x but also on k (or k p ), and even on the polytropic index 
7 for the polytropic case. Thus, we consider the isothermal and polytropic cases separately, 
and obtain the following empirical fits. 

For the isothermal gas, the parameters are fitted to the following polynomials of \i and 

k: 

2 

lo = Uijfi 1 ^ , (37) 

i,j=0 
2 

i,j=0 

s = const., (39) 
2 

t = ^tijiStf. (40) 

i,j=0 

For the polytropic gas, we fix 7 to 7 = 1.15 following Komatsu & Seljak (2001), and fit the 
parameters to the following polynomials of ji and k: 

3 

u = ^W ijf jK? p , (41) 
i,j=0 

3 

C = Yl ^/'''V- ( 42 ) 

i,j=0 

s = const., (43) 
3 

t = Y '<■'„■ (44) 

i,j=0 
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We determine the best-fit values of Uij, Qj, and Uj using the same method as in §2. 
First we calculate E of equation (28) numerically, approximating the integration Ldz' 
to Ldz' where z c = 20i?o which is roughly twice the virial radius of galaxy clusters 
(Makino, Sasaki & Suto 1998). Then we compare the numerical data points of E with the 
fitting formula (36), and determine the best-fit values for each point using the Levenberg- 
Marquardt method. Finally we model the free parameters as polynomials of [i and k (or k p ), 
and determine the best-fit polynomial coefficients. 

The best-fit constant values of s is determined to be unity for both the isothermal and 
polytropic cases just like the 3D model. Figures 3 and 4 plot the best-fit parameters as 
functions of ji for three different values of k (or k p ). It is clear from these figures that the 
polynomial fitting works quite well for all cases. Tables 2 and 3 list the best-fit coefficients 
(ij,Uj,u!ij for the cases of bolometric X-ray and SZ observations, respectively. 

Figure 5 illustrates the degree of the accuracy of our fitting formulae. The filled squares 
represent the numerical results while the solid lines represent our fitting models with the 
best-fit polynomials of (, t, and uj. We have found that our fits reproduce the numerical 
results within the fractional error of 20% for all cases in a range of < A < 10. 

4. DARK HALO RECONSTRUCTION 

In the triaxial dark halo model that we adopt here, the halo reconstruction from the 
X-ray and SZ cluster observations is basically to determine the shapes (e&, e c ) and the ori- 
entations (6, 4>) of the dark halos from the observed surface brightness profiles. 

In §3, we construct a parameterized model for the X-ray and SZ surface brightness 
profiles of galaxy clusters. The model is expressed as a function of the rescaled major axis 
length of the gas isophotes and characterized by three parameters. It turns out that the 
rescaled major axis length and the three parameters depend on the eccentricities and the 
orientation angles of the underlying dark halos, as well as the gas constant. Therefore, we 
find a way to link the observed 2D surface brightness profiles of galaxy clusters to the 3D 
structures of the underlying dark halo and the gas constants of the intra-cluster gas as well. 
In other words, one may expect to determine the values of (e^, e c , 8, 0) and (n, k p ) by fitting 
our 2D parameterized models (eqs. [36] - [41]) to the X-ray and SZ cluster data. 

To demonstrate how our modeling of cluster profiles can be used in the determination of 
the shapes and the orientations of dark matter halos, we apply the following reconstruction 
algorithm to the numerically projected profiles. We first numerically compute £ directly 
using equations (3)- (28). Then we construct a pixeled map of surface brightness in N x N 
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grids (N = 32 in the present example) corresponding to the linear scale of — 10i? < x',y' < 
IORq since 10i? is roughly equal to the cluster virial radius (Jing & Suto 2002). We create 
several realizations of X-ray and SZ profiles for both isothermal and polytropic cases using 
various different values of e&, e c , 9, 0. 

Our reconstruction algorithm proceeds as follows; 

• At each pixel point, say, (x'^y'j), build the model £ th (a^,^) using equation (36). The 
model is characterized by five free parameters e^, e c , 9, 0, k (or k p ) at each point. 

• Fit S th to the observed surface brightness density profiles S obs , and calculate the x 2 '- 



where denotes the observational error at each pixel point (x'^y'j). Those system- 
atic/random errors should depend on specific observation conditions and thus are not 
easy to estimate a priori. Thus in the following numerical tests, we simply assume that 
Oij is unity, independent of i and j. 

• Determine the best-fit values of e&, e c , 9, 0, k (or k p ) through the x 2 minimization. 

The numerical testing results have revealed that the above algorithm works quite well 
in reconstructing the halo eccentricities within the fractional error of 20% as long as the halo 
eccentricities are not so small (e&, e c > 0.3), which can be understood considering that our 
reconstruction algorithm strongly relies on the non-spherical signature, and thus is expected 
to fail for the case of almost spherical halos with low eccentricities (e 6 , e c < 0.3). For the low 
eccentricity cases, the perturbation results of the Paper I may be more useful. Since the low- 
eccentricity case is not the main interest of the present paper, we do not investigate those 
cases here. The algorithm also works in reconstructing the orientation angles but suffer 
relatively large fractional errors. We show the examples of our numerical reconstruction 
results in Figures 6 and 7 for the following cases of the halo eccentricities and orientation 



angles: (e b ,e c ) = (0.4,0.7), (0.5,0.6), and (0.6,0.8); (0,0) = (15, -75), (45, -45), (75,-15) 



in units of degree. 

Figure 6 plots the fractional error of the reconstructed eccentricities versus their input 
values for three different cases of the orientation angles. The left three panels are for e^, 
and the right three panels for e c . The top two panels correspond to the case of 9 = 15° and 
= —75°, the middle two panels to the case of 9 — 45° and = —45°, and the bottom three 
panels to the case of 9 = 75° and = —15°. For almost all cases, the fractional errors fall 
within 20%. 




N 



(45) 



-13- 



Figure 7 plots the fractional error of the reconstructed orientation angles versus the 
original orientation angles for three different cases of the eccentricities. The left three panels 
are for 9, and the right three panels for 0. The top two panels correspond to the case of 
e?, = 0.4 and e c = 0.7, the middle two panels to the case of e& = 0.5 and e c = 0.6, and the 
bottom three panels to the case of e& = 0.6 and e c = 0.8. For most cases, the fractional 
errors fall within 20%. However, for a few cases, the fractional errors are larger than 50%, 
indicating that the angle reconstruction is not always stable compared with the eccentricity 
reconstruction. 



5. DISCUSSIONS AND CONCLUSIONS 

It was regarded as a difficult task to reconstruct the 3D structures of dark matter ha- 
los from the projected 2D surface-brightness maps of X-ray and/or SZ clusters given the 
parameter-degeneracy caused by the projection process itself. To break the degeneracy, pre- 
vious approaches had to rely on such restrictive assumptions as the cluster axial symmetry 
about the line-of-sight and some priors like weak lensing informations. Here we have devel- 
oped, for the first time, a theoretical framework within which the 3D halo reconstruction is 
possible in principle without such restrictive assumptions and additional priors. 

We derived the density profiles of the intra-cluster gas from the 1st principles, assuming 
that the intra-cluster gas is either isothermal or polytropic, and in hydrostatic equilibrium 
within the gravitational potential of concentric and coaxial triaxial dark matter halos. In 
our theoretical modeling of the intra-cluster gas, the density profiles depend explicitly on 
the intrinsic properties of the underlying halos (the two halo eccentricities, e& and e c ). For 
the case of highly aspherical halos, however, we have found that the gas density profiles 
cannot be expressed in closed analytic forms. Therefore, we have attempted to find general 
fitting formulae for the gas density profiles as functions of the halo eccentricities which may 
be applicable even to highly aspherical clusters. We have found empirically a simple scaling 
relation with respect to fx = ejj + e^, and provided a set of fitting formulae for the 3D gas 
density profiles expressed in terms of ji. 

Then, we have incorporated the projection effect into the gas density profiles, and 
provided another set of fitting formulae for 2D X-ray and/or SZ surface brightness density 
profiles of galaxy clusters. Unlike the 3D formulae, the 2D formulae are expressed in terms 
of not only [i but also the orientation angles of the line-of-sight in the halo principal axes. 
In other words, we have found a link of the 2D X-ray and/or SZ observables to the shapes 
and the orientations of the dark matter halos. 
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We have proposed a numerical algorithm based on our model to determine the halo ec- 
centricities and orientation angles from the observed X-ray and SZ surface brightness density 
profiles of galaxy clusters. Ideally we have to apply our algorithm to real observation data, 
but that is beyond the scope of the present paper since it should involve careful treatment 
of the data image analysis together with various observational systematic effects. Thus we 
have decided to test the reconstruction algorithm against simple numerical profiles. We 
found that the algorithm can reconstruct the halo eccentricities fairly accurately if the halo 
eccentricities are greater than 0.3. On the other hand, it turns out that the reconstruction 
of the halo orientation angles are not always accurate, showing large scatters. 

Finally, we conclude that our hydrostatic-equilibrium model for intra-cluster gas in 
triaxial dark halos makes it possible in principle to reconstruct the 3D structures of the dark 
halos from the X-ray and/or SZ cluster maps, as long as the intra-cluster gas and the dark 
halos are well approximated to be in hydrostatic equilibrium and the triaxial ellipsoids with 
a constant axis ratio, respectively. We plan to test our model against real observational data, 
and hope to report the results elsewhere in the future. 
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Table 1: Best-fit polynomial coefficients for the profiles of halo gravitational potentials 





Isothermal 


Polytropic 






P q V 


(3 q 


7] 





1.329 0.426 0.146 


0.780 0.617 


-0.002 


i 


0.118 -0.246 0.036 


0.470 -0.487 


0.133 


2 


-0.010 0.044 0.009 


-0.087 0.136 


-0.013 
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Table 2: Best-fit polynomial coefficients for the profiles of the surface brightness in the 
isothermal case. 







X-ray 






SZ 




(/i re) 
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UJ 


C 




UJ 





3.228 


-1.127 


-0.121 


-3.140 


0.447 


0.014 


1 


-0.206 


0.368 


0.044 


1.133 


-0.087 


-0.005 


2 


0.008 


0.026 


-0.002 


-0.058 


0.018 


0.001 


1 


-0.850 


3.749 


-0.010 
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-0.298 


0.056 


1 1 


1.354 


-1.314 


-0.032 


-1.786 


0.101 


-0.015 


1 2 


-0.126 


0.037 


0.003 


0.166 


-0.027 


0.000 


2 


-12.23 


-2.184 


0.158 


0.297 


-0.003 


-0.042 


2 1 


2.253 


0.840 


-0.025 


0.342 


-0.015 


0.014 


2 2 


-0.082 


-0.044 


0.000 


-0.064 


0.011 


-0.001 
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Table 3: Best-fit polynomial coefficients for the profiles of the surface brightness in the 
polytropic case with 7 = 1.15. 







X-ray 






SZ 




{fi K p ) 
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UJ 
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UJ 





2.261 
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-2.689 
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-0.483 
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-0.049 
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-19.12 


0.424 


1 1 
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-112.8 


-2.301 


-0.685 


77.61 
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1 2 
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0.845 


1 3 
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-14.10 


32.19 


0.145 
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-61.22 


-67.95 


0.836 


36.73 
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-1.279 


2 1 


197.2 
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-1.902 


-125.8 


-147.3 


4.857 


2 2 


-199.1 


-331.6 


0.485 


129.1 


201.9 


-5.537 


2 3 


63.66 


139.9 


0.440 


-39.35 


-79.83 


1.742 


3 


56.23 


38.90 


-0.874 


-28.90 


-18.29 


0.806 


3 1 


-201.0 


-157.9 


2.863 


107.2 


78.73 


-3.186 


3 2 


227.3 
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-2.703 


-123.9 


-110.4 


3.918 


3 3 


-82.35 


-86.20 
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44.58 
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Fig. 1. — Parameters describing our empirical model of the halo potential profile (eq.[19]) 
as a function of the halo eccentricities (p = e\ + e^). The filled circles represent the best- 
fit values to the numerical results at each //, and the solid curves show the corresponding 
polynomial fitting curves. The intra-cluster gas is assumed to be isothermal (Left) and 
polytropic (Right). 
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Fig. 2. — Profiles of gravitational potentials of dark halos. Solid and dashed lines indicate the 
results of the direct numerical integration and our empirical fitting model, respectively. The 
intra-cluster gas is assumed to be isothermal (Left) and polytropic (Right). The eccentricities 
of the underlying dark halos are Top: e& = 0.1 and e c = 0.3, Middle: e& = 0.1 and e c = 0.7, 
Bottom: e& = 0.7 and e c = 0.9. 
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Fig. 3. — Parameters describing our empirical model of the surface brightness profile (eq.[36]) 
as a function of /z for the isothermal case. The filled circles represent the best-fit values to 
the numerical results at each /i, and the solid curves show the corresponding polynomial 
fitting curves. Left: X-ray. Right: SZ. 



-23- 




Fig. 4. — Same as Figure 3 but for the polytropic case. 
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Fig. 5. — Profiles of the cluster surface brightness. The filled circles are computed numeri- 
cally, while our fitting models are shown in solid lines. From left to right panels, we adopt 
the underlying halo eccentricities of (e&, e c ) = (0.1, 0.3), (0.5,0.5), (0.0, 0.7) and (0.6, 0.8). 
From top to bottom panels, we show the cases of isothermal X-ray (k = 8), isothermal SZ 
(ft = 8), polytropic X-ray (7 = 1.15 and K p = 1), and polytropic SZ (7 = 1.15 and n p = 1) 
observations. 
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Fig. 6. — Fractional errors of the reconstructed halo eccentricities. Filled squares, open 
squares, filled circles, and open circles correspond to the cases of isothermal X-ray, isothermal 
SZ, polytropic X-ray, and polytropic SZ, respectively. We adopt the orientation angles of the 
line-of-sight with respect to the cluster in the halo principal frame as (6,4>) = (15°, —75°), 
(45°, —45°), and (75°, —15°) from top to bottom panels. 
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Fig. 7. — Fractional errors of the reconstructed orientation angles. Filled squares, open 
squares, filled circles, and open circles correspond to the cases of isothermal X-ray, isothermal 
SZ, polytropic X-ray, and polytropic SZ, respectively. We adopt the halo eccentricities of 
(e&, e c ) = = (0.4, 0.7), (0.5, 0.6), and (0.6, 0.8) from top to bottom panels. 



